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ABSTRACT 

We compile optical and UV spectra of a sample of 'typical' short period black hole X-ray transients in 
outburst. We also survey determinations of interstellar extinction and distance in order to deredden 
spectra and compare absolute fluxes. Hence we perform a comparative study of the broad-band 
spectral energy distributions (SED). We find that given such a homogeneous sample of typical sources, 
the optical SEDs form a relatively uniform set, all exhibiting quasi-power-law spectra with i^^ oc i/", 
where 0.5 ^ a ^ 1.5 (steeper than the canonical v^^"^ disk spectrum). All become flatter in the UV, 
although there is more diversity here. The SEDs studied can be broadly divided into two optical- 
UV spectral states. The UV-hard spectra, e.g. A 0620-00 and GS 1124-684, continue to rise in the 
far-UV. The UV-soft spectra, e.g. GRO J0422-I-32, drop off. XTE J1859-h226 evolved from UV-soft 
to UV-hard as it decayed indicating that this effect is a real difference, not a dereddening artifact. 
All of the spectra can be fitted by a generalized black body disk model with two forms of heating, 
resulting in the two states. The UV-soft state is consistent with a disk illuminated by a central point 
source, with irradiative heating dominating over viscous. The UV-hard state is well described by a 
viscously heated disk, although this requires very high mass flow rates in the case of GS 1124-684. 
Alternatively, a UV-hard spectrum can be produced if the disk is illuminated by a vertically extended 
X-ray source such as a central scattering corona or jet. Since scattering is assumed by some numerical 
simulations, it is worth emphasizing that where illumination comes from (non-local) scattering high 
above the disk, we generically expect a steeper radial dependence of X-ray heating {F oc R~^) than 
is usually assumed; it is this steep dependence that leads to the UV-hard spectrum. 
Subject headings: accretion, accretion disks - binaries: close ~ stars: individual: V616 Mon, GU Mus, 
V2293 Oph, V518 Per, MM Vel, V406 Vul 



1. INTRODUCTION 

Black hole X-ray transients (BHXRTs) , also known as 
X-ray novae and soft X-ray transients, are low-mass X- 
ray binaries (LMXBs) which undergo dramatic X-ray, op- 
tical, and radio outbursts, separate d by years, decades or 
even longer periods of quiescence (jTanaka fc Shibazakil 
[1996 ChcrcDashchuk 2000). The class is diverse, with or- 
bital periods of 4hrs to a few days to a month. Outburst 
behavior is similarly varied. Among these, however, are a 
sub-group which show relatively consistent behavior and 
these are thought of as the 'typical' members of the class. 
They are often compared to the prototypical BHXRT, 
A 0620-00 (see iKuulkerslll998l for a review). In such a 
typical out burst, the X-ray emission is usually (but not 
always, see lBrocksonp. Bandvonadhvav. fc Fendeill2004D 
dominated by thermal emission from the hot inner ac- 
cretion disk, and UV and optical emission is thought to 
be produced by reprocessing of X-rays by the outer disk. 

A number of attempts have been made to char- 
acterize the spectral energy distributions (SED s) of 
these BHXRTs in outburst. iCheng et all 1)19921) and 
^vncs ct al. (2002) attempted to use simple disk spec- 
tral models, and other authors have used a variety of 
black bodies and power-laws. These mostly only con- 
sider a single source at a time, use a variety of extinc- 
tion curves and reddening determination methods, and 
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do not use comparable spectral models. Consequently it 
is difficult to compare the data on different sources and 
attempt to draw broad conclusions about the population 
as a whole. The goal of this paper is to bring together 
all of the useful data on a suitable homogeneous sample 
of sources, with recalibrations where appropriate for the 
satellite data. It is then possible to deredden the spectra 
consistently and draw direct empirical conclusions about 
the observed SEDs. We can also go on to compare the 
complete sample with simple models in a consistent way. 
Finally we can hope that this compilation will provide a 
reference set as more sophisticated disk spectral models 
are developed. 

Studies of these broad-band SEDs are important as 
they can inform us about the temperature structure 
within the binary. They are especially useful for the 
typical BHXRT systems considered, as these mostly 
have very small mass ratios, hence the emission from 
reprocessing on the companion star is weak, and the 
systems so far discovered have intermediate or low 
inclinations, so disk self obscuration and limb darkening 
should not be large effects. In these cases the com- 
parison between observed SED and disk temperature 
structure should be relatively direct. One must of course 
assume an emergent spectrum for the disk. At present, 
black bodies are probably as good an approximation 
as any for optically thick disks in X-ray binaries, but 
we can anticipate models including a more rigorous 
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TABLE 1 

Sources included in the short period sample considered. 



Object 


Outburst peak 


-^orb 


Adopted 


Dust map 


Distance'' 




(approx. ) 


(hrs) 


E{B - y)=' 


E{B - V)"" 


(kpc) 


A 0620-00 = V616 Mon 


1975 Aug 13 (MJD 42637) 


7.8 


0.35 ±0.02 


0.49 


1.16 ±0.11 


GS 1124-684 = GU Mus 


1991 Jan 15 (MJD 48271) 


10.4 


0.30 ±0.06 


0.37 


5.89 ±0.26 


GRO J0422+32 = V518 Per 


1992 Aug 15 (MJD 48848) 


5.1 


0.35 ±0.10 


0.36 


2.49 ± 0.30 


GRS 1009-45 = MM Vel 


1993 Sep 15 (MJD 49245) 


6.9 


0.20 ±0.05 


0.19 


3.82 ±0.27 


GRS 1716-249 = V2293 Oph 


1993 Sep 30 (MJD 49260) 


<14.7'' 


0.90 ±0.20 


1.08 


(0.24 ±0.04) 


XTE J1859+226 = V406 Vul 


1999 Oct 16 (MJD 51467) 


^ 9.1 


0.58 ±0.12 


0.49 


(8 ±3) 



"Sources for reddening and distance estimates are discussed in the text. 
''Promlgchlcecl. Fi nkbeiner. fc DavisI 1199^^ 

"^GRS 1716-249 does not have a confirmed period, but possible superhumps seen in outburst with a 14.7 hr 
period suggest that the orbital period is similar to this. 



treatment of the disk atmosphere in the future. The 
observed SED essentiaUy can tell us the distribution 
of emitting area as a function of temperature. Given 
assumptions such as axisymmetry and a monotonic 
dependence of temperature on radius then one can map 
this area-temperature relation into the radial temper- 
ature distribution, but the more generalized statement 
is useful in considering more complex geometries which 
could be imagined, such as disks which are warped 
or exhibit spiral waves. In this more general case, 
given a three-dimensional disk geometry, a model for 
the irradiation geometry (e.g. an isotropically emitting 
point source, or perhaps a flat central disk), and some 
method to calculate emergent spectra (at its simplest, 
a black body), one can test specific models against 
the observed spectra, and falsify some combinations. 
There is certainly more than enough uncertainty about 
the temperature structure of irradiated disks to war- 
rant seeking further observational constraints. For 
example, theoretica l work suggests that disks s hould 
be _self-shielding fMever & Mever-Hofmeister Il982t 
:hman. Mincshige. & Wheeler 1990; Dubu s et al.l 
1999I) . and hence that irradiation requires scatterin g 



of X-rays or warping of the disk (jPubus et alJll99^ . 
These two cases may produce very different temper- 
ature structures. An understanding of irradiation 
geometry has immediate application in providing input 
to numerical models of BHXRT outbursts, as these 
models have to make some assumptions about the 
strength, and radial depend e nce, o f X-ray heating (e. 
Kim. Wheeler, k M ineshigej' '1 9991: lEsin et all 1200 



Dubus. Hameurv. fc Lasota 200l UTruss et alJl200^ .' 

We begin in Section |21 by defining our sample of 
BHXRTs. In Section Owe review and evaluate redden- 
ing determinations for these sources, in particular in- 
corporating a reanalysis of 2175 A bump data for some 
sources. This section also discusses the importance of 
the extinction curve adopted. Section 0] discusses corre- 
sponding distance estimates. We then proceed to review 
the database of spectra available based on satellite ob- 
servations (mainly UV) in Section |31 and ground-based 
observations (optical) in Section Having collated our 
sample spectra, we then compare the optical spectra in 
Section [7| and the broader-band optical-UV spectra in 
Section |S| We next take an interlude to evaluate the ef- 
fect of reddening uncertainties upon the spectra in Sec- 



tion IHl before comparing the spectra with disk models in 
Section Uni In Section ITTl we compare our results with 
other studies and with theoretical work. Finally in Sec- 
tion ll2l we summarize our conclusions and consider where 
progress can most effectively be made in the future. 

2. SOURCE SELECTION 

We restrict our sample of sources in order to obtain 
a relatively homogeneous set of spectra. Our selection 
criteria are as follows. 

Firstly we consider only sources which are believed to 
harbor a black hole. In most cases this selection is based 
on a dynamical mass estimate, but where this is not avail- 
able, more subjective criteria are used. Any source ex- 
hibiting type I X-ray bursts is obviously excluded. 

Secondly we exclude those sources with an orbital pe- 
riod greater than one day. Most long period objects ex- 
hibit unusual outbursts, several have companions lumi- 
nous enough to contribute even during outburst and the 
large disks expected in these systems may well behave 
rather differently to those in the short period systems. 

Finally we only include high luminosity systems show- 
ing a 'canonical' lightcurve, i.e. one with a fast rise 
and quasi-exponential decay. Specifically we exclude 
XTE Jl 118+480; this seems appropriate as this source 
exhibited an unusually faint outburst and the optical 
light probably includes significant non-disk emission , 
likely of a synchrotron origin (e.g. iHvnes et al.l 120001) . 
For the sources included, only spectra from the exponen- 
tial decay phase of the lightcurve are used, i.e. we exclude 
mini-outburst spectra. 

An additional practical criterion is of course that there 
must exist flux-calibrated spectra. Based on these cri- 
teria we obtain a primary sample of sources with rela- 
tively complete optical and UV coverage and a secondary 
sample with only optical spectra. Our primary sam- 
ple includes A 0620-00, GS 1124-684, GRO J0422+32, 
and XTE J1859+226. Our secondary sample adds 
GRS 1009-45, and GRS 1716-249. We initially included 
GS 2000+25 in the sec ondary sample , as a fluxed spec- 
trum was presented bv lGharles et alJ l)1991[) and it oth- 
erwise meets our selection criteria. When we came to 
estimate the reddening, however, it was clear that it was 
too large, and too uncertain, for a meaningful derivation 
of the intrinsic spectral energy distribution, so we have 
not attempted an analysis of this spectrum. For all of 
the other targets selected E{B — V)<1 and is relatively 
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well determined. All of these objects also lie at least 5° 
from the Galactic plane. 

For the sources considered, the optical and UV emis- 
sion should be dominated by the disk, although the IR 
may in som e cases include a synchrotr on component (e.g. 
iBrocksoDD et al. ,2002: Corbel fc Fender..2002,) . so we do 
not include IR coverage. 

3. REDDENING ESTIMATES 
3.1. Background 

Since most BHXRTs lie at significant distances, they 
are typically rather reddened objects, hence a reliable 
correction for interstellar reddening is crucial if we are to 
accurately reconstruct the spectral energy distribution. 
This correction depends upon two things: the determi- 
nation of the reddening value, typically parameterized in 
E{B — V), and an extinction curve specifying the amount 
of extinction as a function of wavelength. For a review 
of the problem sec Fitzpatr ick ^1999). 

A number of methods are available to estimate the 
reddening, but they vary considerably in their reliability. 
The most direct measures are not necessarily the best; 
these include the optical color (e.g. B — V) or emission 
line ratios. These do directly measure the differential ex- 
tinction as a function of wavelength, but their interpre- 
tation is model dependent. Less direct methods can be 
more independent of assumptions about the target, but 
instead involve assumed properties for the line-of-sight 
material, for example the gas to dust ratio, the abun- 
dance of Nal, or the strength of the 2 175 A absorption 
feature relative to the broader structure in the extinction 
curve. 

Several approaches also exist for choosing an extinc- 
tion curve, as discussed by Hvnes et al. (2002). As we 
typically have no independent information on the prop- 
erties of the in terstellar rn e dium toward a BHXRT, we 
will adopt the iFitzpatrickl l)1999|) Ry = 3.1 extinction 
curve as this should be the 'best' current Galactic av- 
erage curve. In analyzing our results we must then re- 
member that our uncertainties depend not only on the 
uncertain amount of extinction but the uncertain shape 
of the extinction curve. 

3.2. Methods 
3.2.1. Interstellar spectral features 

Our preferred method to estimate the amount of ex- 
tinction is to measure the strength of the 2175 A inter- 
stellar absorption feature. In fitting the feature we must 
assume some underlying spectral shape, either fixed or 
with some free parameters, and then adjust E{B — V) to 
fit the data. It should be borne in mind that the strength 
of the feature is one of the variable characteristics of the 
UV extinction curve; hence the correlation between the 
2175 A feature and E{B — V) is not perfect, and the in- 
trinsic variance in their relative values introduces an un- 
avoidable 20 p ercent error in reddenings determined in 
this way ( Fitzpatrickl(T999j) . This uncertainty can only 
be avoided if the relative strength of the feature for the 
line of sight is known independently from nearby stars, 
whi ch is only the c ase for one of our sources, A 0620- 
00 l|Wu et al.l 11983(1 . In spite of these limitations, the 
method is generally more precise than the alternatives. 

The equivalent width of the NaD interstellar dou- 



blet can in principle be used to estimate the reddening, 
and calibrations exist for this purpose. The problem is 
that severa l calibrations exist and these do not agree. 
iBarbon et al. (199(|) derive E(B - V) = 0.25 x EW 
from o bservati ons of superno vae in external galaxies. 
IDella Valle fc Duerbeckl i|199g ) on the other hand ob- 
tained E{B -V) = .61 X EW - 0.08 by analyzing the 
data of iCohei] l|197d f) on southern supergiants. Further- 
more, as noted by [Munari & Zwittcr (1997), the NaD 
lines are not in general sensitive to E{B — V) > 0.5 due 
to saturation, so this method can underestimate the red- 
dening; whether it does or not depends on the detailed 
substructure of the lines which is not usually resolved. 
Some caution should therefore be exercised in applying 
this method. 

In addition to the NaD lines, relations between the 
equivalent widths of diff use interstellar bands and E{B — 
V) are commonly used. iHerbid (|1975|) performed an ex- 
tensive study and tabulated relations for 17 interstellar 
bands in the 4400-6850 A range. The correlation be- 
tween different lines was found to be very good but there 
were significant regional effects. 

3.2.2. X-ray absorption 

Another approach taken is to adopt scaling ratios be- 
tween interstellar gas and dust densities, usually ex- 
pressed in a regression relation between iVu and Ay or 
E{B — V). This method is plagued with uncertainty, 
both due to intrinsic variance in this ratio, and due to 
different measures of A'h; it m ay be measured directly 
thoug h UV Ly a absorption ('e.g. lSavaee fc Jenkinsll972 
Bohlin . Savage . & Drake 197^ or indirectly through soft 
X-ray absorption which also includes metallic absorption 
(e.g.' .Gorenstein 1975: Predehl fc Schmitt,1995) : for the 
latter it is necessary to assume abundances. We wi ll 
adopt the relation of IBohlin. Savage, fc Drak^ ((1978(1 : 
— 5.8 X 10^-^ cm~^ x E{B — V), as this relation is rela- 
tively direct and conveniently falls intermediate between 
the other estimates referred to. These authors note that 
individual sources typically scatter by 30 percent around 
this line, and adopting this uncertainty encompasses the 
other average relationships as well. An additional con- 
cern for measurements of Nh is the possible presence of 
local absorption intrinsic to the source. The local absorb- 
ing material will probably be too hot to include a dusty 
component, so will not contribute to the optical/UV ab- 
sorption. 

3.2.3. Outburst colors and line ratios 

Some authors assume a typical intrinsic (B—V) color in 
outburst, and hence derive a direct estimate of E{B — V). 
For our purposes, however, this assumption is extremely 
dangerous; to use such estimates without independent 
corroboration would lead to a circular logic. We would 
be assuming that the intrinsic spectra of all BHXRTs are 
the same to determine their reddening and hence derive 
their intrinsic spectrum. 

A direct measure of the reddening of the spectrum can 
alternatively be obtained from the Balmer decrement, by 
assuming line ratios app ropriate for Case B recombina- 
tion fe.g. lOsterb rock 1989). This approach is obviously 
dependent upon the assumption that Case B line ratios 
are appropriate, and so is also not a very secure method. 
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3.2.4. Extinction maps 

A final technique, of more use for extragalactic 
work, is to use maps of Galactic dust, for example 
ISchlegel. Finkbeiner. fc Davis! l)1998j) . The problem with 
applying this technique to Galactic objects is that it 
only yields the total amount of reddening along that 
line of sight, some of which may be behind the tar- 
get. This technique is thus mainly useful as an upper 
limit, although some sources such as GRO J0422-I-32 
are believed to lie above the dusty component of the 
Galactic disk. An additional concern is that these 
maps do not resolve fine extinction structure and so are 
only r eliable for sources above or below the Galactic 
plane; ISchlegel. Finkbeiner. fc Davis! l)1998j) recommend 
only trusting |6| > 5°. Fortunately, the nature of our 
sample is such that highly reddened objects in the plane 
have been excluded, and all of our targets happen to lie 
outside of this band. We therefore apply this technique, 
remembering that it may include background extinction. 
In most cases the reddening estimated in this way is com- 
parable to, or somewhat larger than, other methods so 
it does seem to be of use, at least for high latitude X-ray 
binaries. 

3.3. Objects 

3.3.1. A 0620-00 

IWu et all llT97^ derived E{B - V) = 0.39 ± 0.02 for 
A 0620-00 from fitti ng ANS UV data through 2175 A fea- 
ture. Subsequently, iWu et all l|1983f l used spectra from 
25 nearby early type stars to derive a line-of-sight ex- 
tinction curve for A 0620-00. It was found that the 
2175 A feature was rather stronger along this line-of- 
sight than average, and hence they derived a revised 
lower estimate of E{B — V) — 0.35. The uncertainty in 
this estimate will be similar to the earlier measurement; 
since they used a line-of-sight extinction curve there is 
no need to account for a 20 percent e rror due to uncer- 
tainty in the intrinsic 2175 A strength. lOke fc Greenstei'nl 
l|1977ft estimated E{B - F) - 0.44, with a large uncer- 
tainty, from the NaD lines. They note, however, that 
the DI Bs seem unus u ally s trong for the reddening de- 
rived. iWhelan et"al] l|1977l) collated the DIB measure- 
ments and confirmed this conclusion. They emphasize, 
however, that the strength of DIBs as a function of 
E{B — V) does vary with Galactic longitude, and that 
A 0620-00 lies within 5° of the region where the DIBs 
are relatively strongest. Using appropriate relations for 
this longitude they derive E{B - V) ^ 0.25 ± 0.05 or 
0.38 ±0.07, depending on the relation used. The redden- 
ing derived from optical interstellar features is therefore 
broadly consistent with that obtained from the UV. Both 
are lower than the value of 0.49 derived from dust maps 
((gchlcgcl. Finkbeiner. fc Davis 1998) but A 0620-00 is 
the closest of our objects and background extinction is 
expected. We adopt the ANS value of 0.35 ± 0.02 for the 
remainder of this work. 

3.3.2. GS 1124-684 

GS 1124-684 was among the best studied 
in the UV with a series of lUE observations 
(jShrader fc Gonzalez-Riestrall993|) early in the outburst 
and a later, higher signal-to-noise HST one l)Cheng et al.l 



I1992D . These per mit fitting of the 2175 A fea t ure to 
derive E{B - V). iShrader fc Gonzalez-Riestral l|1993l^ 
derive E(B - V) = 0.30 ± 0.05 using this method and 
IGheng et all llT99l obtain E{B - V) ^ 0.287 ± 0.004 
from considerably higher signal-to-noise data. Both use 
the Seaton (1919.) extinctio n law. Repea.ting t he analysis 
of the HST data using the iFitzpatricS l)1999(l extinction 
curve, and a more realistic estimate of the accuracy of 
the method, we derive E{B - V) = 0.30 ± 0.06 from 
these data in Section 15.31 consistent with both of the 
earlier estimates. From i n terste llar absorption lines, 
iDella Valle. .Tarvis. fc Wesii l)1991() obtain two values : 
0.3 from applying the calibration of lBarbon et alJ |1990') 
to the NaD lines (with a measured equivalent width of 
I.4A) and 0.35 using the cahbration of iHerbid l|l97,5f) 
on the 5778, 5780, 5797A absorption bands. These are 
consistent with the UV determination. The hydrogen 
column density has been estimated from Ginqa data 
alone as about 1.6 x 10^^ cm^^ (Ebisawa et al. 1994), 
and from Ginqa ajuA ROSAT spectra at 2.2 x 10^^ cm^^ 
(jGreiner et al. Il994|) . These imply reddenings of 
E{B — F) ~ 0.3 and 0.4 respectively, which are just 
about consistent with optical estimates given the 
uncertainty in the relation between the two. For the re- 
mainder of this work we will adopt our revised estimate 
based on the 2175 A feature: E{B -V) ^ 0.30 ± 0.06. 
As for A 0620-00, this is slightly lower than was derived 
from dust maps, 0.37. 

3.3.3. GRO JO422+32 

The 2175 A feature in GRO J0422+32 has been mea- 
sured by Shrader et al. ( 1994) using lUE data. They 
obtain E{B - U) = 0.40 ± 0.06. We have repeated the 
analysis using a more recent extinction curve and ob- 
tain a range of values of 0.29-0.48, depending on the 
spectral range modeled and the underlying spectrum 
which is assumed (see Section |^2J . Many authors have 
also attempted to derive a red dening from interstellar 
hues. IHarlaftis fc Charles! (IT991 s uggest E{B-V ) = 0.2 
based on 4428 A and 6613 A DIBs. IShrader et all l|1994f) 
also obtain E{B - V) = .2 ± 0.1 from the 5780A 
DIB. Finally iCallanan et"an ((199 5) measure equivalent 
widths of 1.2 ± 0.2 A for the NaD lines and 0.5 ± 0.1 A 
for the 5780/5788/5797A group. Using_ the relations 
of Barbon et al. (1990) and^Hcrbig ( 197al) they ded uce 
E{B -V) = 0.3 ± 0.1. IChevalier fc Ilovaiskvl lll99l ex- 
amine 12 stars in the field surrounding GRO J0422-I-32 
and find that their color-color diagram needs to be dered- 
dened by E{B - V) = 0.40 ± 0.07 in order to match 
the unreddened main sequence. This method assumes 
that all these stars, and GRO J0422-I-32, are subject to 
the same reddening; it would be appropriate if all are 
above the main plane of extinction. Even if some are 
not, but GRO J0422+32 is (as seems likely given the 
estimated distance for GRO J0422-I-32), th en this value 
becomes a lower limit. IShrader et alJ l)1997j) fit absorbed 
power-laws to ASGA spectra from the 1993 August mini- 
outburst to obtain iVn — 1.6 x lO^^cm"^, correspond- 
ing to E{B -V) = .28 ± 0.08. Finally, and most 
recently, iGelino fc Harr ison (2003) derive E{B — V) = 
0.24 ± 0.03 from fitting the quiescent spectral energy dis- 
tribution, but this value is somewhat sensitive to assump- 
tions about the amount and color of disk light. There 
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is clearly little close consensus on the reddening value, 
and even the method which is usually most precise, fit- 
ting the 2175 A feature, yields equivocal results. Conse- 
quently, based on the spread in measured values we adopt 
E{B -V) ^ 0.35 ± 0.10. This choice is rather arbitrary, 
but is representative of the scatter between the derived 
estimates. It is completely consistent with the estimate 
of 0.36 from dust maps. 

3.3.4. GRS 1009-45 

ine11a Valle eTaTI iTmfll present several arguments for 
a low reddening in GRS 1009-45, E{B - F) - 0.2. 
From the strength of NaD absorption they derive an 
estimate of E{B ~ V) — 0.18, with an upper limit of 
E{B -V) = 0.23. Assuming that the intrinsic {B - V) 
color is similar to other LMXBs implies a consistent value 
of E{B — V) ^ 0.23 (although see caution above). We 
follow these authors in adopting E{B-V) = 0.20 ±0.05. 
The column density is not well constrained, but proba- 
bly < 1.2x1021 cm-2 ||Kubota et al. 1998), which implies 
E{B — V) < 0.27, consistent with the optical estimates. 
Finally dust maps predict E{B — V) < 0.19, consistent 
with other estimates. 

3.3.5. GRS 1716-249 

The reddening of GRS 1716 249 has been estimated by 
[Delia Vallc. Mirabel. & Rodriguez (1994) as E{B-V) = 
0.9 ± 0.2 based on a number of arguments giving roughly 
consistent values. The strength of the NaD lines implies 
E{B — V) > 0.75 (see the note on saturation of these hnes 
above). The strength of DIBs suggests E{B — V)= 0.6- 
0.9. Arguments based on the outburst (B — V) color, and 
the Balmer decrement give somewhat higher values, but 
are model dependent as discussed above. The interstel- 
lar hydrogen column density has been estimated at A^h = 
4 X 10^1 cm-2 (Tanaka 1993), implying E{B -V) = 0.5- 
0.9, which is consistent with the optical estimates. Fi- 
nally, dust maps imply E{B — V) 1.08. In summary, 
we see no compelling reason to amend the estimate of 
iDella Valle. Mirabel, fc Rodriguez! l)1994j) . leaving this as 
the largest and most uncertain reddening in our sample. 

3.3.6. XTE J 1859+226 

iHvnes et all pOCil estimated E{B - V) = 0.58 ± 
0.12 using th e well defined 2175 A interstellar feature. 
iWagner et all (199?) obtained E{B - V) = 0.8 ± 0.4 
from the 5780 A DIB. Estimates of the hvdrogen column 
density included A^h = 0.8 x 10^^ cni"^ ijdal Fiume et all 
[1999I from SAX data) and A^h = 1-1 x cm-^ 
lMarkwardil l200H RXTE). This column density corre- 
sponds to E{B — V) ~ 0.9-2.5, allowing for the uncer- 
tainty and variance in the gas-dust ratio. Even with this 
uncertainty, and that in the 2175 A determination, there 
is then a discrepancy. In contrast, dust maps imply an 
even lower value, E{B ~ V) < 0.49. The 2175 A value 
thus represents a reasonable middle-ground, so we will 
use t his estimate , and the derived dereddened spectra 
from IHvnes etaTI l|200^ . 

4. DISTANCE ESTIMATES 

The determination of distances to BHXRTs is even 
more problematic than their reddenings, since there are 
no direct measures of distance. Ideally, we would seek 



some kind of standard candle to determine the distance. 
For neutron star LMXBs, a standard candle does exist 
in the form of radius expansion type-I X-ray bursts. No 
such analog exists for black holes, and any attempt to use 
the outburst luminosity, as done by some authors, runs 
into a circular logic problem again, that we would be ar- 
tificially forcing uniformity on our sample before compar- 
ing them. The only reasonably robust method that has 
been applied to our sample is to use the quiescent com- 
panion star. System parameters provide estimates of the 
companion's temperature and radius, and hence abso- 
lute magnitude, albeit with a large uncertainty in many 
cases. Given the observed magnitude, estimated extinc- 
tion and the fraction of quiescent light due to the disk 
(the veiling) we can then derive a distance. The errors 
are likely to be large, but are representative of the real 
uncertainty in knowing the distance, and other methods 
are likely to be even worse. The best method to assess 
the errors is to perform a simple simulation, constructing 
a population of model binaries consistent to within the 
observational uncertainties and derive a distance for each 
one, and hence a mean distance and standard deviation. 
This method has recently been applied to a number of 
our sources using IR data. This is preferable because 
the IR derived distance is less sensitive to uncertainty 
in extinction or companion temperature than the opti- 
cal is. The best estimate is for A 0620-00 for which 
iGelino. Harrison fc Or^ i200lD derive 1.16 ± 0.11 kpc. 
For GS 1124-684 lGelino. Harrison, fc McNama ra (200]J|) 
quote ~ 5.1 kpc with no uncertainty, but this estimate 
was subsequently refined with simulations as above to 
5.89 ± 0.26 kpc (jG elino 2002). For GRO J0422-I-32 
iGehno fc HarrisonT 42003.) estimate 2.49 ± 0.30 k pc and 
finally for GRS 1 009-45 iGelino fc HarrisonI l|200^ derive 
3.82 ± 0.27 kpc. Of course, these are not the only esti- 
mates for these sources. However we believe they are the 
most thorough and reflect the most up to date parameter 
estimates. These estimates are broadly consistent with 
previous values. Some concerns do still exist for these ob- 
jects in that these estimates, and their quoted uncertain- 
ties, depend on the limits placed upon the veiling. The 
works cited generally assume rather small veilings in the 
IR, and while th is assumption is widely made, its v alidity 
is questionable l|Hvnes. Robinson, fc BitneJ^l2004^ . The 
most immediate effect of increasing the veiling is to re- 
duce the apparent brightness of the companion star, im- 
plying a larger distance. It is therefore quite possible that 
the distances are underestimated s omewhat. For more 
detailed discussion of these issues see lJonker fc Nelemand 
iDOi). 

For the other two of our selected sources, the infor- 
mation necessary to use the quiescent companion as a 
standard candle is simply not all available yet, so in these 
cases we make use of other indicators from the literature. 

The distance to G RS 171 6-249 is the poor est deter- 
mined in our sample. IDella V allc. Mi rabel, fc Rodrigue3 
l) 19941) estimate 2 kpc from relations between Na D equiv- 
alent width and distance. Since the NaD lines could 
well be saturated for the relatively high extinction of 
this object, this distance is really only a lower limit, 
even assuming that the mean relation is valid along 
this line of sight. Based on the peak optical bright- 
ness, these authors estimate an upper limit of 2.8 kpc, 
but as noted above we should be cautious about using 
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such measurements for our purposes. In this case, how- 
eve r, we have httle other information t o go o n and fol- 
fow iDella Valle. Mirabel. &: Rodriguez! 1)19941) in adopt- 
ing 2.4±0.4kpc for GRS 1716-249. The true uncertainty 
may well be larger, however, and the distance could be 
significantly greater. 

The distance to XTE J1859+226 is also problem- 
atic. iZurita et al.l l)2002|) estimated d = llkpc from 
both the brightness of the outburst and that of the 
quiescent companion. The latter estimate was, how- 
ever, based largely upon guesses for a number of bi- 
nary parameters incl uding the orbital period and c om- 
panion spectral type. iFilippenk o fc Chornocld ({20011 ) re- 
port spectrosco py indic ating an o rbital period longer 
than assumed bv'^Zurit a'et al.l (|2002(1 and an earlier spec- 
tral type. Both effects will tend to make the qui- 
escent companion more luminous and hence push the 
source to a larger distance. We have done a simula- 
tion similar to those discussed above to test these ef- 
fects, and estimate that the quiescent brightness implies 
a distance of ~ 21 ifc 5 kpc as suming the constraints of 
iFilippcnko fc Chornockl 1)200 ID . This value depends on 
some uncertain or totally unknown parameters, however, 
such as the inclination and mass ratio. The disk veil- 
ing is also ignored and a significant disk contribution 
would push the distance even higher. The faintness in 
quiescence therefore suggests a very large distance, large 
enough to conflict with other estimates. Outburst X- 
ray observations place a likely constraint of 5-13 kpc, 
based on a combination of spectral and timing informa- 
tion (^arkwardt 2001). Fits to the optical -UV spectral 
energ y distribution suggested 4.6-8.0 kpc l)Hvnes et al.l 
120021) ■ assuming a plausible range of possible disk radii 
and inclinations. In view of the extreme uncertainty 
about the system parameters, in particular the compan- 
ion star's effective temperature, we believe these outburst 
estimates are currently more reliable, but that we have 
no really trustworthy distance for this object. We adopt 
a working estimate of 8±3kpc, which encompasses most 
of the parameter space proposed thus far, but note that 
a larger distance remains a possibility. 

5. SATELLITE SPECTRA 

5.1. ANS 

A 0620-00 was observed by the Ultraviolet Experiment 
onboard the Astronomical Netherlands Satellite (ANS) 
around day 47 after X-ray peak (Wu et al. 19 76); a sec- 
ond observation was made around day 225 ijWu et al.l 
ri983l) which we do not use as it was obtained b eyond the 
maiii exponential decay phase of the outburst. l)Wu et al.l 
11976)) tabulate the data from the earlier observation. 
From these we constructed an average using only the six 
observations when the broad 1550 A bandpass was used 
to minimize the impact of the possible C IV emission line. 

5.2. lUE 

lUE spectra w ere obtained for two o f our sources: 
GS 1124-684 ( S hrader & Gonzalez-R iestra 1993) and 
GRO J0422-h32 jShrader et all 11994)) . Both sources 
were observed with the low-resolution short- and long- 
wavelength prime cameras (SWP and LWP). We ob- 
tained all of these data fr om the archive, using the 
NEWSIPS reprocessed results l)Nichols fc Linskvlll996ft . 



For GS 1124-684, good coverage was obtained through- 
out the outburst, although the later data are of poorer 
quality. We assembled three representative average spec- 
tra. For the early outburst, we combined spectra from 
1991 January 17-24; no obvious spectral differences are 
visible between these. We produced a second epoch 
spectrum from February 6-28. This corresponds to 
some way down the decline, but before the secondary 
maximum in the X-ray lightcurve. For the third sam- 
ple we use a single epoch, March 27-28, within the 
secondary maximum, wh ich was seen in the UV data 
l)Shrader fc Gonzalez-Ries tra 1993). At later epochs, the 
lUE data are of poor quality, and the HST spectrum (see 
next section) is of more use. We did not attempt to re- 
derive a reddening for GS 1124-684 from the lUE data 
as the HST spectrum of this source is of much higher 
quality. 

For GRO J0422-t-32, the coverage was less complete; 
the first two months were observed regularly, but there 
was only one more observation about a year later. We 
assembled two average early spectra from 1992 August 
16-26 and 1992 September 2-27 respectively. The last 
visit lies beyond the main exponential dec a.v phase of the 
outbu rst. Since the reddening derived bv lShrader et alJ 
()1994)) from the 2175 A feature was inconsistent with 
other estimates, we combined these two spectra and re- 
analyzed them to a ttempt to resolve the discrepancy. 
IShrader et aI\E99^ denveA E( B -V) = 0.40±0.06, but 
did so by assuming the entire lUE spectrum was a power- 
law. Consequently there is a danger that the derived 
reddening may be biased by a fai lure o f the model. As 
for XTE J1859-I-226 <)Hvnes et aL))200l(l . we tried using 
a range of different wavelength ranges and both a power- 
law and irradiated disk model in order to test for such 
sensitivities. We do indeed find that a wide range of red- 
dening values are derived, 0.29-0.48, dependent mainly 
on th e choice of fitt i ng ran ge. As can be seen from the 
fit bv lShrader et al.l l)1994)) . a power-law model does fail 
both at short and long wavelengths. We find that using 
the 1650-3120 A range adopted for XTE J1859-I-226, sig- 
nificant systematic residuals are seen at the longer wave- 
lengths. Fits with this region (A > 2500 A) included pro- 
duce large values of E{B — V) = 0.37 — 0.48, whereas if 
this region is excluded a range of 0.29-0.38 is found. The 
range above 2500 A contains little if any information on 
the 2175 A feature, and the discrepancy with other es- 
timates is also largely resolved if it is excluded, so we 
choose to use fits without this region. 

5.3. HST 

HST spectra from the main part of the outburst 
are availab le for two of our sources: GS 1124-684 
(iCheng et al. 1992) and XTE J1859-t-226 l)Hvnes et alJ 
l200ir 

The data on GS 1124-684 were obtained with the 
Faint Object Spectrograph (FOS) fairly late in the out- 
burst, but still in the exponential decay phase. We 
obtained these data from the ST-ECF archive. They 
were rec alibrated by the ar chive using the P0A_CALF0S 
pipeline l)Alexov et al.ll200lj) and the final reference files. 
As for other source s, we reanalyzed the 2175 A feature. 
ICheng et all (|1992) found a good fit using a power-law 
model and we find the same; the irradiated disk model 
which worked better for XTE J1859-I-226 provides a very 
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TABLE 2 

Satellite spectra included in the sample. 



Object 


Date 


Outburst 


Facility 


Wavelength 


Plotted'' 


Reference'' 






Day 




range (A) 






A 0620-00 


1975 Sep 28-30 


46-48 


ANS 


1550-3300 


Y 


1 


rm 1 1 9zi— Rj^A 


1QQ1 Tan 17—91 


2—6 


lUE 






2 




1991 Jan 24 


9 


lUE 


1150-3300 


A 


2 




1991 Jan 29 


14 


lUE 


1150-3300 


N 


2 




1991 Feb 6 


22 


lUE 


1150-3300 


A 


2 




1991 Feb 12 


28 


lUE 


1150-3300 


A 


2 




1991 Feb 18 


34 


lUE 


1150-3300 


A 


2 




1991 Feb 25 


41 


lUE 


1150-3300 


A 


2 




1991 Feb 28 


44 


lUE 


1150-3300 


A 


2 




1991 Mar 15 


59 


lUE 


1150-3300 


N 


2 




1991 Mar 27-28 


71-72 


lUE 


1150-3300 


N 


2 




1991 Apr 22 


97 


lUE 


1150-3300 


N 


2 




1991 May 15 


120 


HST/FOS 


1580-4810 


Y 


3 


GRO J0422+32 


1992 Aug 16 


2 


lUE 


1150-3300 


A 


4 




1992 Aug 20 


6 


lUE 


1150-3300 


A 


4 




1992 Aug 26 


10 


lUE 


1150-3300 


A 


4 




1992 Sep 2 


19 


lUE 


1150-3300 


A 


4 




1992 Sep 10 


27 


lUE 


1150-3300 


A 


4 




1992 Sep 20 


37 


lUE 


1150-3300 


A 


4 




1992 Sep 27 


44 


lUE 


1150-3300 


A 


4 


XTE J1859+226 


1999 Oct 18 


2 


HST/STIS 


1150-10 200 


Y 


5 




1999 Oct 27 


11 


HST/STIS 


1150-10 200 


Y 


5 




1999 Nov 6 


21 


HST/STIS 


1150-10 200 


Y 


5 




2000 Feb 8 


115 


HST/STIS 


1150-10 200 


Y 


5 




2000 Mar 5 


141 


HST/STIS 


1150-10 200 


Y 


5 



'^Spectra shown in Fig.Qor Fig.U] Y indicates a spectrum plotted individually, A means the spectrum 
is indued in an average and N that the spectrum is omitted. 

''References: 1. Wu et al. (1976), 2. Shrader & Gonzalez-Riestra (1993), 3. Cheng et al. (1992), 4. 
Shrader et al. (1994), 5. Hynes et al. (2002) 



poor fit for GS 1124-684. A reddened power-law fit to the 
1574-3500A range (chosen to exclude any weak Balmer 
jump contribution) gave E{B — V) — 0.30, a — 0.52. 
The re ddening estimate is not significantly different to 
that of iCheng et alJ l)1992|) . but the derived spectrum is 
somewhat bluer. This is d ue to the dif f erent extinction 
curve used: we use that of Fitzoatrick (1999), whereas 
ICheng et all l)1992|) used the .Seaton 1,1979) curve. The 
differe nce between these curves is shown bv iFitzpatrickl 
(jl999(l . and the different spectral slopes derived is repre- 
sentative of the real uncertainty in the process (see Sec- 
tion]^ for more discussion of this point.) We adopt this 
value, with 20 percent uncertainty, as a reddening mea- 
sure for GS 1124-684: E{B - V)'^ 0.30 ± 0.06. We also 
repeated the fit with the full dataset, 1574-4816 A, and 
with all the wavelength ranges used for XTE J1859-I-226. 
In every case, we derive E{B — V) consistent with the 
range E{B -V) = 0.30 ± 0.06. 

A full analysis of data obtained with the Space Tele- 
scope Imaging Spectrograph (STIS) of XTE J1859-H226 
along the lines here has already been presented by 
iHvnes etTaLl l|2002f) . we simply include selected results 
of that work for comparison with the other sources pre- 
sented here. These spectra span the full 1150-10 300 A 
range at several epochs through the outburst. These are 
the only data on this source included. 

6. GROUND-BASED OPTICAL SPECTRA 
6.1. A 0620-00 



The highest qual ity calibrated spec t ra for A 0620-00 
were presented by lOke fc Greensteinl l)1977j) who show 
spectrophotometry from several epochs, of which two, 
around 21 and 118 days after outburst peak respectively, 
were obtained during the exponential decay phase of the 
outburst. The original data are no longer available, so 
these two spectra were digitized from the paper and in- 
cluded in our sample. The first ob servation occurred 
close to the ANS UV observations. iGull et all l|1976n 
present two low-resolution SEDs of the optical. These are 
of poorer quality than those of Oke & Grecnstcin ( 1977), 
and are only shown with an arbitrary fiux scale, so we 
do not include these. 

6.2. GS 1124-684 

A large sample o f optical spectr a were pre- 
sented by [Pc lla Val le. Jarvis. & WestI l|199lD and 
iDella Valle. Masetti. fc Bianchinil l|1998t) . These span 
days 0-124 days after outburst peak The spectrum 
included in IDella Valle. Jarvis. fc West! (|l991) was dig- 
itized; those from the later paper were kindly supplied 
by the authors. For the latter, the reddening correction 
pre-applied was removed before dereddening them as 
described above for consistency. For further analysis, 
we merged the spectra from 1991 Jan 18-23. We also 
rejected the 1991 Mar 25 spectrum which has limited 
wavelength coverage, and a very different slope to all 
others considered, even 1991 Mar 21. Some portions at 
the extreme near-IR and near-UV ends of some spectra 
show abrupt and non-repeatable transitions and these 
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TABLE 3 

Ground-based spectra included in the sample. 



Object 


Date 


Outburst 


Facility 


Wavelength 


Plotted'' 


Reference'' 






day 




range used (A) 






A 0620-00 


1975 Sep 2-4 


20-22 


Palomar, Hale 5 m 


3200-10 350 


Y 


1 




1975 Dec 8-11 


117-120 


Palomar, Hale 5 m 


3200-10 350 


Y 


1 




1 QQ1 Ton 1 ^ 

iyyj_ J (ill ±u 


n 
u 


TT^O NTT 




Y 


2 




1991 Jan 18 


3 


ESQ, 1.5m 


3800-8600 


A 


3 




1991 Jan 19 


4 


ESQ, 1.5m 


3800-7400 


A 


3 




1991 Jan 20 


5 


ESQ, 1.5m 


6300-8600 


A 


3 




1991 Jan 23 


8 


ESQ, NTT 


4200-6200 


A 


3 




1991 Feb 15 


31 


ESQ, 1.5m 


3850-7450 


Y 


3 




1991 Feb 21 


37 


ESQ, 1.5m 


3850-9850 


N 


3 




1991 Mar 6 


50 


ESQ, NTT 


4000-8150 


Y 


3 




1991 Mar 21 


65 


ESQ, 2.2m 


4050-9050 


N 


3 




1991 Mar 25 


69 


ESQ, NTT 


6100-8100 


N 


3 




1991 May 5 


110 


ESQ, 1.5m 


4000-7100 


N 


3 




1 QQ1 A/Tnv 1 Q 


124 


hjOW, o.uiii 




Y 


3 


GRO J0422+32 


1992 Aug 17 


3 


Lowell Perkins, 1.8 m 


3800-8500 


Y 


4 




1992 Aug 18 


4 


Lowell Perkins, 1.8 m 


4850-6800 


N 


4 




1992 Sep 25 


42 


Lowell Perkins, 1.8 m 


4850-6800 


N 


4 




1992 Oct 10 


57 


Lowell Perkins, 1.8 m 


4850-6800 


Y 


4 




1992 Nov 4 


82 


Lowell Perkins, 1.8 m 


4850-6800 


N 


4 




1992 Dec 15 


123 


Lowell Perkins, 1.8 m 


4850-6800 


Y 


4 




1993 Jan 21 


160 


Lowell Perkins, 1.8 m 


4850-6800 


N 


4 




1993 Feb 11 


181 


Lowell Perkins, 1.8 m 


4850-6800 


Y 


4 


GRS 1009-45 


1993 Nov 18 


64 


ESQ, 2.2m 


3600-9050 


Y 


5 


GRS 1716-249 


1993 Oct 5-8 


5-8 


ESQ, 1.5m & 2.2m 2 


4600-6700 


Y 


6 



^Spectra shown in Fig. Y indicates a spectrum plotted individually, A means the spectrum is indued in an 
average and N that the spectrum is omitted. 

''References: 1. 'O ke fc GreensteinI fwrHl. 2. IDella Valle. Jarvis. &: WesH <T99ll). 3. 

Ipella Valle, Masetti, & Bianchini (1993), IShrader etaU <1994l) . 5^ IDella Valle et al.l 119971) . 6. 

Ibella Valle. Mirabel, Rodriguez, a994l^ 



were also excluded. 

6.3. GRO J0422+32 

IShrader et alJ l|1994|) present spectra from 1992 August 
16 to September 25. We include an absolutely fluxed 
spectrum from 1992 August 17, from this dataset, which 
was provided by the author. The remaining epochs were 
digitized from a figure also provided by the author as the 
original data were not available. 

6.4. GRS 1009^45 

GRS 1009-45 suffered both from a rel atively short ex- 
ponential phase of outburst, ^ 150 days l|Bailvn & Orosd 
Il995'l . and from a late discovery. The only spec- 
tr um from the exponentia l decay phase was presented 
bv IDella Valle et alJ l(T997|) . spanning 3500-9500 A. This 
spectrum was digitized, placed on an absolute flux scale 
(assuming 1 unit = 10~^^ ergcm"^ s~^; Delia Valle, pri- 
vate communication), and corrected to our preferred 
dereddening solution. 

6.5. GRS 1716-249 

One flux calibrated spectrum from the ex- 

ponential decay phase wa s presented by 

IDella Valle. Mirabel, fc Rodrisrued l|1994) . This spec- 
trum was digitized and corrected to our preferred 
reddening value and extinction curve as for other 
sources. 

7. THE OPTICAL SPECTRA 



Our combined sample of optical spectra are shown in 
Fig. n For these purposes, satellite spectra have been 
truncated at log:^ = 14.98 Hz, the approximate highest 
value recorded for A 0620-00 from the ground. We have 
excluded some spectra for which the flux calibration ap- 
pears incorrect; when there are large distortions which 
are not present at any other epoch; and others to avoid 
duplication of two very similar spectra. 

It is clear that after applying the source selection cri- 
teria described above, and subject to these additional 
restrictions, the spectra define a relatively homogeneous 
sample. All are blue, with F^, increasing toward the UV. 
All have an approximately power-law form. Subject to 
the large uncertainties in distance described above, all 
represent comparable optical luminosities. 

To perform a more quantitative comparison, we must 
parameterize and fit the spectra. At this stage we con- 
sider only very simple characterizations; we will calculate 
model disk spectral energy distributions in Section 1101 
but UV coverage is required to adequately constrain 
such models. The two simple characterizations widely 
used in the literature are a power-law or a hot black 
body; the latter has an approximately power-law form 
on the Rayleigh-Jeans tail, and so can also provide an 
adequate fit to the data. Given that a single temper- 
ature black body is likely a very poor description of an 
accretion disk, the black body fit has limited physical sig- 
nificance, but is useful for comparison with other work. 
The temperatures that we derive are somewhat lower 
than other authors have obtained in this way. In the 
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Fig. 1. — Optical spectral energy distributions for all of our sources. Open and filled symbols are used to distinguish alternate epochs 
and have no other significance. Numbers indicate the approximate number of days after X-ray outburst peak when each spectrum was 
obtained. The bow ties indicate the 1 it uncertainty in spectral slope due to uncertain dereddening of the spectra. The point with an 
error bar to the right of each panel indicates the uncertainty in absolute calibration for each object due to uncertain distance and absolute 
extinction. Solid lines are power-law fits. 
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TABLE 4 

Results of fits to the optical spectra. 



Object 


Outburst 


Power-law^ 


T~> 1 „ „1 "U „ J 

Black body 




day 


a 




1 (K) 


A 0620-00 


20-22 


0.94 ± 0.01 


±0.07 


14 700±500 




117-120 


0.81 ±0.02 


±0.07 


12 900±500 


GS 1124-684 





n -t- n nfi 
u.oo in u.uu 


^u.zz 


14100±800 




3-8 


1 07 -1- n OR 


±0 21 


17100±500 




31 


J. . U -L ^ U. U'i 


±0.22 


17 000±400 




37 


0.99 ± 0.04 


±0.20 


15 300±600 




50 


1.20 ±0.03 


±0.22 


19 800±800 




65 


0.60 ±0.03 


±0.21 


11800±300 




liU 


i.4y it U.U4 


±U.zU 


oi oUUitvSUUU 




124 


J-.tJii _1_ IJ.VJrt 


±0.20 


20100±1000 




4 


1.07 ± 0.03 


±0.36 


1 7 700+1 ono 




5 


1.05 ±0.05 


±0.39 






42 


0.57 ± 0.02 


±0.39 


1 3 100+300 




57 


0.57 ±0.04 


±0.39 


1 3 1 00+500 




82 


0.57 ±0.03 


±0.39 


13 100+400 




123 


0.50 ±0.02 


±0.39 


12 600±300 




160 


0.39 ±0.03 


±0.39 


11800±300 




181 


0.48 ±0.02 


±0.39 


12 400±200 


GRS 1009-45 


64 


1.19 ±0.03 


±0.17 


17 000±700 


GRS 1716-249 


5 


1.07 ±0.07 


±0.78 


17 200±900 


XTE J1859+226 


2 


1.51 ±0.01 


±0.40 


28 000±600 




11 


1.36 ±0.01 


±0.40 


22 500±400 




21 


1.26 ±0.01 


±0.40 


19 900±400 




115 


1.10 ±0.03 


±0.42 


17400±600 




141 


1.11 ±0.05 


±0.42 


17100±900 



''Defined by oc 

''Uncertainty due to uncertain E{B — V). 



case of A 0620-00 we estimate 13 000-15 000 K, compare d 
to 25 000-30 000 K obtained bv lOke fc GreensteinI l(l977|) . 
For both GRS 1716-249 and GRS 1009-45 we estimate a 
temperature of ~ 17 OOP K, compared to p ublished esti- 
mates of 20 000 K JPelTa Valle et al.llT997() and 25 000 K 
l)Della Valle fc Duerbecklll993j) respectively. The differ- 
ences likely reflect the difficulty of reliably constraining 
such high temperatures using optical data alone, and in 
this regime large derived temperature differences can re- 
sult from quite small differences in the data and its treat- 
ment (for example the wavelength range used in fitting.) 
In view of this difficulty, and the limited physical mean- 
ing of a single temperature fit, we ignore the results of 
these fits from hereon; suffice to say that the spectra in- 
dicate a hot continuum, sufficiently hot that optical ob- 
servations have very limited discriminatory power alone. 

Power-law fits have previously been performed for 
some of our sources. They provide a more empiri- 
cal description of the data, but one which more di- 
rect ly relates to its information co ntent. For GS 1124- 
684 iKing. Harrison, fc McNamaral (|T996) found an ini- 
tially blue spectrum, a ~ 0.7 (where F^, cx j/"), which 
softened to a ^ 0.2 after the secondary maximum. 
iCheng et alJ l)1992() also remarked that this object ex- 
hibited a period at the beginning of the outburst when 
the source had a steeper (blue r) spectrum th an later ob- 
served. For GRO J0422+32, IShrader et all llT99l find 
the optical spectral index softened from 0.6 to 0.0 over 
the first three months of the outbur st, then stabilized. 
IKing. Harrison. Sz McNamaral 1)19961) found spectral in- 
dices of between -0.1 and -0.4 for most of the outburst. 



although the evolution is less clear in their data and they 
assumed a relatively low extinction, Ay — 0.72. We fit 
power-laws to our optical spectra and tabulate power- 
law indices in Table 2| These values differ somewhat 
from those quoted above, as absolute values are rather 
sensitive to the assumed E{B — V). Some spectral indices 
have also been plotted as a function of time in Fig. [3 All 
of the values that we derive fall within the approximate 
range +0.5 < a < +1.5, and it should be noted that 
all are steeper (bluer) than the canonical a = 1/3 of a 
steady-state viscously heated disk which is sometimes as- 
sumed. Some variations within this range may be due to 
incorrect dereddening, but the full range cannot be ex- 
plained in this way, and the variations between repeated 
observations of a single source indicate that there are 
real differences. There is no single 'correct' value, and 
so the utility of assuming a standard spectral slope (or 
photometric color) in estimating reddening values is lim- 
ited. The observed variation corresponds to an intrin- 
sic {B - V) color of 0.00 to -0.24. This range is smaU 
enough that it does have some value for an approximate 
reddening estimate for sources with relatively large red- 
dening, but one should still be cautious, as we have only 
obtained this degree of uniformity by a posteriori selec- 
tion of sources by outburst behavior and orbital period. 
For example, the unusual source XTE J1118+480 has 
a somewhat flatter optical spectrum (Hv nes et alJl200fH) 
and the long-perio d system GR O J 1655-40 is much red- 
der in the optical l|Hvnes et alJ lT"9981 

For A 0620-00 there are too few points to draw 
any convincing conclusions about evolution of the op- 
tical spectrum, other than that little change appears 
to have occurred between the two epochs. This con- 
clusion is b orne out on longer times c ales b y the pho- 
tometry of iLlovd. Noble, fc PenstonI l|1977j) who find 
approximately constant B — V values around +0.25 
throughout the outburst, corresponding to a dereddened 
color of —0.1, comparable to our spectral shapes. For 
GS 1124-684 there are more data, but these exhibit 
rather large scatter. The spectra may be hardening, 
but it is hard to say with confidence, and the ghoto- 
met ric analyses of [K ing. Harrison^ fc McNamaral (Il99l 
and lDella Valle etaL (1997.) suggest the opposite behav- 
ior. GRO J0422+32 and XTE J1859+226 show more 
coherent trends, but this may reflect better sampling 
and/or data quality rather than an intrinsic difference. 
The latter in particular, being based solely on satellite 
data, should be robust. Both sources clearly become 
softer optically during the outburst. There appears to 
be rapid evolution during the first month, then slower 
softening (at least for GRO J042 2+32) later. T h is evo - 
lution has already been noted by 'S hrader et alJ (|1994), 
although all of their power-law indices are softer than 
ou rs due to a lower assumed redd e ning. The photometry 
of IKing. Harrison, fc McNamaral l)1996t) may also indi- 
cate an initial rapid softening. This behavior suggests 
that the optically bright regions are cooling during the 
outburst. We will defer more detailed discussion until we 
consider the combined optical and UV spectra of these 
sources (Section^)). Although the evolution is very sim- 
ilar for GRO J0422+32 and XTE J1859+226, the actual 
values are rather different. This difference may be real, 
or may reflect an incorrect dereddening of the spectra. 
The dotted hues plotted in Fig. El are offset by 0.45. At 
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Fig. 2. — Evolution of the optical power-law index (_F^ <x u"). 
The dashed line corresponds to the canonical u^^^ value for a 
steady-state viscously heated disk; all spectra appear steeper than 
this. The dotted line is shows an exponential decay fitting the 
data for XTE J1859+226, and is intended only to guide the eye 
in comparing these data with those for GRO J0422+32, for which 
the line has been offset downward by 0.45. 



late times GRO J0422+32 is about a further 0.25 below 
the plotted line. Uncertainty in reddening introduces 
an uncertainty in the spectral slope of ^ 0.4 for both 
sources, however, so most of the difference in the optical 
slopes could be due to an incorrect E{B — V). We will 
reexamine the difference in the light of the UV spectra 
and models in Section [TUI 

We have invested some time in comparing the optical 
results in isolation as these are often all that is avail- 
able. In addition these results can be compared directly 
with optical photometry which typically has better tem- 
poral sampling. Based on optical results alone, the ob- 
jects appear relatively homogeneous and we cannot be 
confident that apparent differences are not due to in- 
correct dereddening. Optical colors do give some useful 
information about spectral evolution independent of red- 
dening uncertainty. They do not allow us to distinguish 
between different models robustly however; for that we 
must extend the data into the UV. Based on the analysis 
so far, GRO J0422+32 and XTE J1859-f 226 exhibit def- 
inite and similar softening during the outburst, whereas 
A 0620-00 and GS 1124-684 exhibit no clear trends, al- 
though very little data are available for the former, and 



the latter exhibit large scatter. 

8. OPTICAL AND UV SPECTRA 

Our combined sample of broad-band spectra is shown 
in Fig. 13 We have selected data for which both optical 
and UV coverage is available. Where possible we have 
used contemporaneous data, but in some cases, particu- 
larly for A 0620-00, the time difference is significant. 

It can be seen that the power-law fits to the optical 
data do not adequately extrapolate for any object. All 
UV spectra flatten with respect to the optical. The de- 
gree of flattening varies, with A 0620-00 and GS 1124- 
684 continuing to rise more slowly in the UV, while 
GRO J0422-h32 and XTE J1859-h226 actually turn over 
and begin to fall. One might question whether this dif- 
ference is real or reflects incorrect dereddening, but since 
we see both forms in XTE J1859-I-226, it is likely a real 
variation within the sample. We can characterize this 
distinction by dividing the spectra into two optical/UV 
states. We refer to spectra in which the far-UV spectrum 
continues to rise as UV-hard, and ones in which it turns 
over and declines as UV-soft. 

The presence of such a turnover is natural and is a 
consequence of the cutoff in the disk temperature distri- 
bution at the disk edge (or at a cooling front). We will 
discuss the interpretation of these spectra in Section [TUI 
First, however, we examine more carefully how the spec- 
tra are affected by uncertain dereddening. 

9. THE IMPACT OF UNCERTAIN DEREDDENING 

9.1. Optical 

The effect of uncertain dereddening is easiest to quan- 
tify in the optical, as optical extinction curves are less 
diverse than those in the UV, and the variation of extinc- 
tion with wavelength is approximately monotonic. Un- 
certainty in E{B — V) will introduce an uncertainty in the 
optical spectral slope, whereas an uncertain Ay will af- 
fect the implied luminosity level. Note that these effects 
are semi-independent since Ry = Ay/E{B — V) may not 
always have the Galactic average value of 3.1. To quan- 
tify the effect on the spectral slope, for each source we 
have shown in Fig. ^ a bow-tie indicating the power- law 
slopes derived if E{B — V) is allowed to vary within the 
uncertainties estimated. We have also included a corre- 
sponding estimate of the change in power-law index, Aa, 
in Table H 

9.2. UV 

The effect on the UV spectrum is more complex, as 
the UV extinction curve exhibits several features, and 
the relative strength of these can vary. The effect of un- 
certain E{B — V), assuming a fixed extinction curve, can 
be quantified in the same way as for the optical data, 
although the form is not simply a variation in power-law 
slope. This effect is shown by the curved bow-ties in 
Fig. 131 Note that there is no increase in the uncertainty 
around the 2175 A bump as this region is explicitly flat- 
tened out in dereddening the spectra; a larger E{B — V) 
would only be consistent with the data if the bump were 
unusually weak and vice versa. 

Fitzpatrick (1999) has presented a convenient param- 
eterization of the variations in the extinction curve, 
and also presents a recipe for quantifying the ef- 
fect of the resultant uncertainty in the derived SED. 
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Fig. 3. — Broad-band spectral energy distributions for the sources which have UV data. Notation is as for Fig.Q The distorted bow-ties 
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iHvnes et al.l l|2002D have already applied this analysis 
for XTE J1859-I-226. We show comparable results for 
GS 1124-684 and GRO J0422-H32 in Fig.EI Note the un- 
certainty in the extinction curve is less of a problem for 
A 0620-00 for which a line-of-sight curve was determined 
in pa rallel with the reddening determination (,Wu ct al.. 

ITosl . 

10. DISK MODELING 

10.1. The basic model 

We have shown that while a single power-law is an ac- 
ceptable description of most optical spectra, it does not 
adequately describe the UV break, which seems to be 



ubiquitous. This break is to be expected. We beheve 
that the spectra included in this sample should reflect 
the spectrum of an accretion disk in a hot state. The 
hot region should terminate at effective temperatures be- 
tween 6 000K and 10 000 K, as regions at this tempera- 
ture should be unstable. If such a cut-off occurs, then 
the optical spectrum will lie on the Rayleigh-Jeans tail 
for most of the disk area hence exhibits a steeper (bluer) 
spectrum than at higher frequencies, as observed. 

We can attempt to quantify the disk temperature dis- 
tribution using the commonly used model for a disk 
heated by both viscosity and irradiation and emitting 
locally as a black-body. This model, its origins, and orig- 
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inal references, are summarized bv IHvnes et al.l lj2002j) . 
and we will not repeat the discussion. An alternative dis- 
cussion is provided by [Chakrabartv ( 199 ^. We us e that 
model exactly as described bv IHvnes et al.l lj2002tl . and 
with the same definitions and terminology. We empha- 
size that the model actually only describes a disk with 
two forms of heating, one of which drops off as (vis- 
cous heating) and one as approximately (irradiation 
heating). The exact model used (based on lVrtilek et al.l 
Il99 ff) assumes a slightly slower drop-off of the latter and 
is adopted for consistency with earlier work and as a sim- 
ple empirical description, r ather than as a th eoretically 
motivated formulation (see lDubus et alJll999|) . In prac- 
tice the exact prescription adopted has little effect on 
the optical and UV spectra, and so any error introduced 
in this way will be less than that due to uncertainties 
in dereddening the spectra. The essential discriminating 
evidence between these two effects is the shape of the far- 
UV spectrum; a UV-hard spectrum which continues to 
rise favors r~'^ heating, whereas a UV-soft spectrum with 
a clear peak indicates heating. Beyond the peak, the 
UV-soft spectrum will tend towards a spectrum or 
steeper, but this form will not be realized within the 
observable spectral range. The component, which 
gives rise to the classical v^^'^ accretion disk spectrum 
ijShakura fc Sunvae v 19731 , is usually assumed to be due 
to viscous heating, but could equally describe irradiation 
which drops off in a similar way with radius. Indeed this 



functional form is commonly adopted in mode ling irra- 
diated disks in protostars fe.e:. IHartmannllT998|) . We will 
return to this alternative interpretation in Section llU.51 
The analy sis we will perform is virtually identical to 
that done bv lHvnes et"al] l|200^ for XTE J1859-f 226 and 
we will directly compare to those results without reana- 
lyzing those data. One important difference, however, is 
that XTE J1859-H226 is the only source for which a full 
UV-optical spectrum was obtained quasi-simultaneously 
with HST. For the other objects in our sample having 
UV coverage, the optical spectra were taken from the 
ground days or even weeks apart. The non-simultaneity 
compromises the reliability of the SEDs obtained, and 
wc will discuss the severity of this problem on a case- 
by-case basis. In general it requires that we introduce 
an additional free parameter representing the difference 
in flux levels between the two epochs, but it should be 
remembered that spectral changes might have occurred 
as well. 

It is worth emphasizing that we explicitly assume there 
are no other sources of flux besides the accretion disk; 
emission from a jet or heating of the companion star 
are neglected. For the sources and wavelength region 
adopted, this is a reasonable assumption, and was a ma- 
jor criteria in their selection. 

10.2. A 0620-00 

A 0620-00 has a high quality optical SED, and a well 
determined UV extinction, so in principle it should be 
one of the best suited objects to our study. It suffers, 
however, from very crude wavelength coverage in the UV 
and, most severely, from the lack of simultaneous UV 
and optical observations. We will compare optical data 
taken around day 21 of the outburst with UV data from 
around day 47. The 26 day interval between them may 
mean that the spectrum has changed significantly. A 
particular concern is that the UV observation happened 
around the time of the ons et of a seconda ry maximum 
in the X-ray lightcurve f e . g. iKuulkersI 1 998|l . The optical 
spectrum taken later, around day 118, is very similar to 
the early one however, suggesting that little evolution 
actually occurred in the optical spectrum, and that the 
comparison is valid. 

Keeping these concerns in mind, we proceed to fit a 
disk model to the SED, allowing for an arbitrary off- 
set between the optical and UV data. This approach is 
equivalent to assuming that the only difference between 
the two epochs was in the flux level of the spectrum, 
rather than the shape. The fits are shown in Fig. |S1 
The best fit model has no irradiative heating and corre- 
sponds to the classic viscously heated steady-state disk. 
For comparison we also show a purely irradiated model. 
This is an extremely poor fit, even with an arbitrary off- 
set between the optical and UV. Improving the UV fit 
would require an increase in the irradiation temperature 
and would steepen the optical spectrum. Of course some 
optical spectral evolution may have occurred. The later 
spectrum is flatter, however, and other sources also tend 
to flatten as the outburst proceeds. Flattening of the 
optical spectrum would further degrade the irradiated 
fit, so it seems unlikely that the irradiated model can be 
accommodated in this way. 

The implied parameters (for the UV epoch) are sensi- 
ble. The outer edge viscous temperature, 7400±400K, is 
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TABLE 5 

Parameters of disk fits to spectral energy distributions. 



Object 


Dates 


Ti 


rr,out (-^) 


vise 


,out (K) 


Normalization 


A 0620-00 


20-48 





0-6 000 


7400 


7 000-7 800 


11.35 


9.65-13.35 


GS 1124-684 


2-9 





0-6 000 


10 000 


8 600-11400 


0.85 


0.61-1.26 


GS 1124-684 


22-44 





0-5 000 


7 800 


7 000-8 600 


0.98 


0.75-1.31 


GS 1124-684 


120-144 


7 000 


5 000-9 000 


8 400 


8 000-8 800 


0.23 


0.19-0.27 


GRO J0422+32 


2-10 


12 600 


11 400-14 200 


5 200 


3 000-5 600 


1.31 


0.98-1.72 


GRO J0422+32 


19-44 


9 400 


8 400-11000 


5 000 


4 600-5 600 


1.98 


1.28-2.78 
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Fig. 5. — Disk model fits to the spectra of A 0620-00. Since the 
optical and UV data are non-simultaneous, the separation has been 
retained and models have been offset to fit. The solid line shows 
the result of an unconstrained fit, which favors a disk heated only 
by viscosity. The dashed line shows a fit constrained to Tvisc = 0, 
i.e. a purely irradiated disk. The latter clearly cannot fit well. To 
improve the UV fit would require a large increase in the irradiation 
temperature and would then steepen the optical spectrum to an 
unacceptable degree. 



comparable to that found in XTE J1859-f 226 late in out- 
burst, and is close to the expe cted m inimum temperature 
of the hot phase of the disk (Lasota 20011). The normal- 
ization can be interpreted as described by iHvnes et al.l 
((2002). Assuming a distanc e d = 1.16 ± 0.11 kpc, and 
inclination i = 40.75 ± 3.0° (jGelino. Harrison. &: Oros3 
I200 1'). we estimate a disk radius rout = (1-8 ± 0.3) x 
10^^ cm, comparable to the expected tidal truncation ra- 
dius of the disk (1.7 x 10^^ cm). The implied disk radius 
would increase to 2.0 x 10^^ cm when the model is nor- 
malized to fit the optical data. 

We can also estimate the mass flow rate required 
to provide t his o ptical luminosity as described by 
IHvnes et all 1)20021) . At the UV epoch, we estimate 
M = (0.9 ± 0.3) X IQ-'^Mf:, yr-^ Assuming an II.OMq 
black hole l)Gehno. Harrison, k Oro's3l200lj) . and 10% 
accretion efficiency, the Eddington limit corresponds to 
M = 2.5 X lO^'^Mgyr-i, so we obtain a mass flow rate 
a few tenths of the Eddington limit. 

Given the uncertainties, the observations are thus con- 
sistent with a lobe-filling, viscously heated disk. Con- 
trary to common assertions, there seems to be no need for 
irradiative heating to provide the energy radiated in the 



optical-UV. This is not to say that irradiation plays no 
role; it may still be crucial in stabilizing the disk against 
a transition to a cool state. The point is that the main 
source of energy emitted is consistent with viscous heat- 
ing rather than irradiation. 

Finally, we note that the long wavelength spectrum 
is flatter than the model predicts. This difference is 
relatively subtle, and may not be of great significance. 
However, a similar l ong-wavelength flat tening is seen 
in XTE J1859-t-226 f Hvnes et all 12002) and may in- 
dicate a weak contribu tion from synchrotron emission 
ijBrocksopD et al.ll2002() . 

10.3. GS 1124-684 

GS 1124-684 has better temporal coverage than 
A 0620-00, and generally optical observations can be 
selected from within the time-period spanned by mul- 
tiple lUE observations. In addition, the HST spectrum 
obtained later in the outburst is of higher quality, albeit 
lacking the shortest wavelengths. Two concerns are man- 
ifest here, however. Firstly, we have already seen that the 
power-law slopes derived from the optical spectra scatter 
significantly. The most likely interpretation of this scat- 
ter would be that the flux calibration is not totally reli- 
able. If this is correct then the spectral fits are somewhat 
compromised. Secondly there is a sharp upturn in the 
far-UV, rather similar to that seen in XTE J 1859+226 
when the latter was corrected with the pure Fitzpatrickl 
ll'999'l extinction curve. Hvnes et al. (2002) interpreted 
the upturn in that case as due to a far-UV rise compo- 
nent o f extinction which is weaker than in the lFitzpatrickl 
l)1999j) extinction curve. The same could be true for 
GS 1124-684, and similarly introduces uncertainty in 
the dereddened far-UV spectrum. We therefore omitted 
these points from the fits. 

Disk model fits are shown in Fig. (HJ performed as 
for A 0620-00. GS 1124-684 is also best fitted by a 
model dominated by viscous heating, exhibiting UV-hard 
spectra at all epochs. A model with purely irradiative 
heating would not fit. The disk outer edge tempera- 
ture is comparable to that in A 0620-00 at late times, 
but higher earlier, dropping from 10 000 ± 2000 K to 
~ 8 000K; a similar decline was found bv iCheng et alJ 
(,1992i). Assuming d = 5.89 ± 0-26 kpc and i = 54 1^°^ 
l|Gehno. HarrisonTfc McNamarall200lt iGelind l20?)l we 
obtain disk radii of 1.4-3.0x10^^ cm, comparable to the 
expected tidal truncation radius of 1.7 x 10^^ cm. The 
scatter inferred in the radius may arise from the uncer- 
tainty in the optical spectral slope, since this uncertainty 
will affect the derived temperature, which in turn will be 
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Fig. 6.— Disk model fits to the spectra of GS 1124-684. For the 
first two epochs optical and UV data do not overlap or line up well 
so they have been fitted separately. The last epoch spans a shorter 
period, and the data overlap, so no ofli'set was included in this case. 
The solid line shows the result of an unconstrained fit, which favors 
a disk heated only by viscosity. The dashed Une for Epoch 3 shows 



a fit constrained to T„: 



■ 0, i.e. a purely irradiated disk. 



inversely correlated with the normalization. These data 
thus also appear approximately consistent with a lobe 
filling disk. The higher temperatures early on may in- 
dicate that the disk is indeed being tidally truncated, 
since the hot-phase of the disk should be able to extend 
to below 10 000 K if it were unbounded. 

The temperatures we estimate are higher than for 
A 0620-00, even though the disk radius is expected to 
be comparable. A higher temperature implies a larger 
mass flow rate (if the spectrum is indeed due to viscous 
heating). Note that this conclusion is largely indepen- 
dent of distance u ncertainti es, provided both disks are 
tidally truncated. lEsin et a l. (2000) indirectly provide 
support for this conclusion. They compare the X-ray 
lightcurves and estimate a flux scaling factor of ^ 8 
for the ratio of A 0620-00 to GS 1124-684. They argue 
that the overall flux level of the two outbursts should be 
the same, assuming that the Eddington scaled luminos- 
ity of the intermediate state should be constant. This 
is only expected, however, for extreme parameters, for 
example a black hole mass in A 0620-00 of just 4.5 
is assumed. If we instead adopt current best estimates 
of the parameters JGelino . Harri son, fc McNamaral2001t 
iGelino. Harrison fc OrosJ200H) . then we instead predict 
a flux scaling of 77 = 56 ± 18, if the Eddington scaled lu- 
minosities were the same. This scaling is uncertain, but 
significantly larger than observed, implying that A 0620- 
00 is underluminous relative to GS 1124-684, and hence 
that th e accretion rate in GS 1124-684 was likely larger. 
lEsin et al. (2000) also note that the ratio of optical to 
X-ray flux is a factor of 1.6 larger in A 0620-00 than in 
GS 1124-684, and attribute this difference to more effi- 
cient irradiation of the outer disk in A 0620-00. There is 
another explanation, however. Both the optical and X- 
ray bandpasses lie at the extreme ends of the spectral en- 
ergy distribution, on the Rayleigh- Jeans and Wien tails 
respectively. If the mass flow rate is higher, the effect 
is to increase the temperature throughout the disk and 



shift the SED out of the optical bandpass and into the 
X-ray one. Such a difference is supported by the hotter 
optical/UV spectra inferred in GS 1124-684. This spec- 
tral shift will naturally decrease the optical to X-ray flux 
ratio without invoking any other changes. For example, 
consider two viscously heated disks around 10 black 
holes, with inner and outer radii of 3 Rsch and 2 x 10^^ cm 
respectively. We adjust the mass flow rate to give outer 
disk temperatures of 8 000K and 12 000 K respectively 
and calculate the ratio of V band to X-ray flux (1.3- 
7.5 keV). The ratio can change by at least a factor of 
two and possibly more, depending on the relative weight- 
ing across the X-ray bandpass and the effect of spectral 
hardening on the X-ray spectrum. This effect is more 
than adequate to explain the difference in the flux ra- 
tios between A 0620-00 and GS 1124-684 without any 
irradiation effect. 

If we assume a 6.95 Mq black hole in GS 1124- 
684, then we estimate mass flow rates of (2.0 ± 1.2) x 
lO-^Moyr-i, (8.5 ± 1.2) x 10-'^Moyr-\ and (1.3 ± 
1.2) X 10~^Mo yr~^ for e pochs 1 -3 respectively (compa- 
rable to that obtained bv lMisral fl999) assuming similar 
parameters). These values are superficially somewhat 
problematic, as the Eddington limit should be reached 
for M = 1.6 X 10~^ M0 yr~^, assuming 10 % accretion ef- 
ficiency. Even with higher efficiencies, the early phase of 
the outburst should thus have involved super-Eddington 
accretion. This is not a fatal problem, however, since 
the Eddington limit strictly only applies to spherical 
accretion and a geometrically thin disk could deliver 
mildly super-Eddington accretion rates. Such a super- 
Eddington model was indeed cons i dered for GS 1124-684 
by IMisral l|1999(l . and lEsin etsl\ (|2000D also assumed a 



peak accretion rate somewhat above the Eddington limit. 
A more pertinent question is whether such a large accre- 
tion rate can be reconciled with the observed thermal 
(dis k) component i n the X-ray spect rum of GS 1124- 
684 ()Ebisawa et al.lll994(l . IMisral lfT999i) did address this 
question in the context of their spectral modeling and 
show that the two are consistent provided that Compton 
scattering modifies t he X-ray spectrum with a spectral 
hardening factor fShi mura fc Takaharalll995h of fc ^ 3. 
Thi s factor is larger than t he value of 1.7 usually adopted, 
but iMerloni. Fabian, fc R oss ( 2 d00(l argued that values 
as large as 3 are plausible. If one accepts this correc- 
tion then both the optical/UV spectrum, and the soft X- 
ray thermal component can be consistently interpreted 
with a purely viscously heated disk. As noted above for 
A 0620-00, irradiation may still play a role in stabilizing 
the disk in the high viscosity state, but it is not required 
to provide the main energy input. We will consider an 
alternative explanation for the optical/UV spectrum in 

Section nnini 

10.4. GRO JO422+32 

For GRO J0422-I-32 we have two useful epochs with 
both lUE and optical data. These data may also suffer 
from similar problems with the far-UV rise component 
of extinction to GS 1124-684 and GRO J0422-H32, but 
in the opposite sense, as there seems to be an enhanced 
drop off in flux at the shortest wavelengths. 

Fits to these data are shown in Fig. [3 Strong irra- 
diative heating is needed in this case, although as the 
outburst declines viscous heating may be becoming im- 
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Fig. 7.— Disk model fits to the spectra of GRO J0422+32. 
Off^sets have been included to allow for non-simultaneity, but only 
a small difference is required. 



portant as in XTE J1859-I-226; the timebase and data 
quality are inadequate to say this with any confidence, 
however. The irradiation temperatures inferred are lower 
than for XTE J1859+226, with ~ 13 000 K early in 
outburst, com pared to - 18 000 K in XTE J1859+226 
l|Hvnes et al.l 120 02'). Lower temperatures are consistent 
with the softer optical spectrum earlier estimated. As 
seen in Fig. 01 this difference does seem to be reflected 
in the lower peak frequency of the SED, so is likely 
real. As for the other sources, we can estimate the 
disk radius r equired; assuming d = 2 49 ± 0.30 kpc and 
z = 45 ± 2° l|Gelino fc Harrisonll200l . we derived disk 
radii of (1.3 ± 0.3) x 10" cm and (1.7 ± 0.5) x 10" cm 
for the first and second epochs respectively. These num- 
bers seem too large, as we estimate a tidal truncation 
radius of ^ 9 x 10^° cm assuming the parameters of 
IGelino fc HarrisonI l)2003f) . This discrepancy could be re- 
solved with a reduction in the source dist ance (which is 
not totally impossible; iWebb et al.l l2000^). with a more 
modest increase in the temperature (which is factored 
into the large uncertainty in this derived radius), or if 
the inclination uncertainty is larger than assumed. 

10.5. Indirect heating of the disk 

We have seen that UV-hard spectra are ubiquitous, 
and can occur even when a BHXRT is at a high optical 
luminosity. It is common, at least in the field of X-ray 
binaries, to interpret such a spectrum as a signature of 
a vis cously heated steady-state disk (e.g. iCheng et al.l 
11992(1 . and thus far we have done the same. This is not 
implausible, particularly in the case of A 0620-00 for 
which the required mass flow rate is not large. 

All this spectral form really indicates, however, is 
the temperature distribution, and that the heating rate 
varies as approximately . An alternative irradiation 
model is used by the young star community which ac- 
tually predicts irradiative heating varying as (e.g. 
Eartmann 1998). This model addresses the case where 
a circumstellar disk is irradiated by the stellar surface. 
The key difference is that since the star is not compact, 
it is not a point source, but an extended one which il- 



luminates the disk from above. Indeed a similar case 
which also produces heating arises from a 'lamp- 
post' model, where the so urce of X-rays are raised above 
the disk fe.g. iDumonTfc CoUin-Souffrin 1990). 

As described in the appendix, one can produce UV- 
hard spectra with an elevated (or vertically extended) 
X-ray source. This behavior is generically expected, ex- 
cept for the specific case of a concave disk with a suffi- 
ciently compact X-ray source, possibly mediated by lo- 
cal scattering by an extended disk atmosphere. It is 
worth emphasizing that where irradiation is attributed 
to X-rays scattered by a central coronal source, one ex- 
pects heating, in contrast to the assumption usu- 
ally made in modeling outburst behavio r that irradia- 
tion drops off approximately as R^^ (e.g. King & Ritted 
1998, Kim. Wheeler, & Mineshigc 1999, Esin ct al. 20qJ 
iDubus. Hameurv. fc Lasotal 120011) . If X -ray scattering 
occurs in a highly vertically extended disk atmosphere, 
however, i?""^ heating could still be appropriate. 

How plausible are these scenarios? X-ray scat- 
tering 'coronae' are certainly believed to be present 
in X-ray binaries. However, it is likely that coro- 
nae are not spherical, but very flattened, and repre- 
sent a vertical extension of the d isk atmosphere (e.g. 
iChurch fc Balucihska-Churchll200 4'). Effective illumina- 
tion of a convex disk, with a vertically extended central 
corona, thus seems difficult to achieve, and hence this 
seems unlikely as a mechanism for producing UV-hard 
spectra. 

Finally, we note that X-rays from a jet might provide 
an alternative 'lamppost' illumination, cither by direct 
synchrotron emission or by scattering. Relativistic beam- 
ing might then be expected to considerably reduce the 
back-illumination of the disk, but large luminosities are 
not required, especially as the X-ray source could then 
be higher above the disk than is plausible for a corona. 

10.6. Heating of a warped disk 

IDubus et alJ l(1999t) suggested that warping of disks 
may facilitate irradiation of disks which would oth- 
erwise be self-shielded, and this could also be ex- 
pect ed to modify the radial dependence of the heat- 
ing. iQgilvie fc DubusI lf2001„) have examined this ques- 
tion quantitatively. They do find strong variations in the 
radially averaged heating, but these are not as simple as 
the additional 1 / R dependence that would be required to 
reproduce the F^, oc v^l'^ spectra, although they do find 
that an average reprocessing efficiency of C = 5 x 10"'^ 
can readily be produced with a warped disk. While only 
a fraction of the disk is illuminated at any time, the ef- 
fect of irradiation on the thermal structure of the disk 
will last for of order the local thermal timescale. This 
should be comparable to the local dynamical timescale 
(i.e. the Keplerian rotation period), so heating effects will 
be azimuthally averaged as material rotates through the 
illuminated patch and hence the whole disk will be af- 
fected (G. Dubus, priv. comm.). Further investigation of 
the spectral energy di stributions expected in the warped 
disk case is needed. IQgilvie fc~ D ubus (2001) do note, 
however, that the short-period BHXRTs we are consid- 
ering here are expected to be stable against radiation 
driven warping, so in the absence of another warping 
mechanism, it remains questionable whether warping can 
be an important factor for these objects. 
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11. DISCUSSION 

In light of the above discussion, we should reconsider 
the origin of the optical emission in BHXRTs in outburst. 
It is sometimes stated that optical emission is always to- 
tally dominated by reprocessed X-rays. This may be true 
in neutron star systems, but there is good reason to ex- 
pect irradiation to be weaker in the black hole case, as 
there is no direct emission from the com pact object sur- 
face ijKing. Kolb. fc Szuszkiewic3ll997(l . This difference 
is invoked by these authors to explain the observations 
that black hole LMXBs are transient, whereas most anal- 
ogous neutron star systems are persistent X-ray sources. 

The studies presented here suggest that the situation 
may be more marginal than usually assumed. If one in- 
terprets a J^^/^ spectrum as a signature of a viscously 
heated disk, then the answer is clearly no, as this spec- 
trum occurs rather often among the population. In the 
case of A 0620-00, at least, if one takes a plausible mass 
flow rate, and plausible disk size then viscous heating 
alone can readily explain the observed optical brightness. 
In the case of XTE J1859-I-226, the SED can initially be 
interpreted as irradiation dominated, but a transition to 
the viscously dominated regime appears to occur suggest- 
ing that the irradiation dominance is rather marginal. 
Of course this conclusion hinges upon one interpretation 
of the SED, but the alternatives considered are relatively 
unpalatable. Irradiation by X-rays from a lamppost (per- 
haps the jet?), or by a vertically extended corona could 
reproduce the observed SED, but requires an excessive 
luminosity for the extended source and/or a rather large 
vertical extent. We thus appear left with viscous heating 
of the outer disk as the most natural explanation of the 
observed optical flux during outburst, at least when UV- 
hard spectra are seen. As noted several times, irradiation 
should still play a role in maintaining the disk in a hot, 
high viscosity state by suppressing the disk instability, 
but it does not appear to dominate the energy budget. 

The point is often made that studies such as this one 
are severely compromised by assumption of local black- 
body spectra. At some level, of course, this is true, 
but in the UV, and especially the optical, the effect is 
likely to be mild. The two most important effects on a 
optical-UV broad-band SED, for the temperatures con- 
sidered, are hydrogen breaks and limb-darkening. Line- 
blanketing by many blended absorption lines could have 
some effect, although the highest quality data suggest a 
rather clean continuum where the effect of intrinsic (non- 
interstellar) absorption lines is very weak or absent, ex- 
cept for Lya (Haswell et al. 200^. The Balmer jump 
falls well within our spectral coverage. While this re- 
sults in strong distortion of the SED of hot stars, such 
large effects are clearly not present in our data, and there 
is little if any evidence for its presence at all. The ob- 
servations are consistent with spectral models of accre- 
tion disks around stellar mass black holes which do pre- 
dict rather weak Balmer jumps (e.g. iHubenv et aLll200H 
Hubeny priv. comm.) A much larger effect is predicted 
at the Lyman break, and has the effect of significantly 
enhancing the far-UV flux longward of the break. This 
effect produces an upturn which is principle is detectable. 
Such an upturn wa s indeed observed in XTE J1859-I-226 
l|Hvnes et al.l2002j) . and may also be present in GS 1124- 
684, as shown earlier. An alternative explanation was ad- 



vanced bv lHvnes et al.l l)2002(l . however, that this upturn 
originates from deviation of the far-UV rise component 
of extinction from the Galactic average. Disentangling 
the two effects is likely impossible without line-of-sight 
extinction curves or a larger sample of objects. Fur- 
thermore, strong Lyman breaks can als o be predicted in 
AGN spectra, but are not observed (see lBlaes et al.ll2'001l 
and references therein). The other deviation from black 
body spectra that could be si gniflcant is limb-darkening 
ijDiaz. Wade, fc Hubenvl IT996 ') . This is wavelength de- 
pendent and stronger in the UV, so tends to soften the 
far-UV flux. It is an inclination dependent effect, how- 
ever, and for the moderate inclinations of the objects 
considered here (40-60°; see Sectional should not be 
dominant. With both Balmer and Lyman breaks and 
limb-darkening, the models of Hub env et al.l l)200lD do 
not deviate from their black-body counterparts by more 
than ±0.1 in logF^ over the optical-UV range, for incli- 
nations of less than 60°. In particular, a UV-hard spec- 
trum remains UV-hard. Consequently black-body spec- 
tra do appear adequate for the analysis presented here, 
and it is likely that the dominant uncertainty is actu- 
ally in the reddening corrections, not the local emergent 
spectra. 

12. CONCLUSIONS 

The optical and UV spectral energy distributions of 
BHXRTs in outburst have the potential to inform us of 
the temperature distribution of the outer disk, and hence 
test models for the heating of these disks. We have col- 
lated all the spectra published to date on short period 
typical systems which we believe are useful in this regard. 
In spite of a 25 year baseline, these data are rather lim- 
ited. HST data are by far the most suitable for this work 
in both the UV and the optical. Ground-based optical 
studies suffer from uncertainties in the flux calibration, 
especially at short wavelengths, and from Telluric ab- 
sorption at long wavelengths. Consequently even from 
the small sample available, only a subset are really reli- 
able. 

Even with high quality data, further complications 
arise from uncertainty in the line-of-sight extinction 
curve. If the curve is known, then the reddening can 
be determined to high-precision using the 2175 A feature. 
Without a priori knowledge of the extinction curve, how- 
ever, the relative strength of the feature is not known and 
the reddening cannot be precisely determined. Further 
difficulties arise from variations in other aspects of the 
UV extinction curve. The far-UV rise seems particu- 
larly variable in our sample and several of the sources 
considered exhibited kinks in the far-UV which could be 
attributed to an incorrect correction of the far-UV rise. 
These further compound the uncertainty in the dered- 
dened spectral energy distributions. 

In spite of these difficulties, we can draw some use- 
ful conclusions. Provided a suitable sample of typical 
short-period BHXRTs are selected, the optical spectra 
form a relatively homogeneous sample. All have quasi- 
power-law shapes, with spectral slopes (in F^, vs. i^) of 
0.5-1.5. Where the data are good enough to discrimi- 
nate, the optical spectra tend to soften as the outburst 
decays, corresponding to the cooling of the outer disk. 
The very homogeneity of the optical spectra means that 
they are of limited value alone in testing disk models, and 
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they cannot discriminate between different forms of disk 
heating. It is also worth remarking that since the optical 
lies on the long wavelength tail of the spectral energy 
distribution the V band flux is not necessarily related in 
a straightforward way to the irradiating luminosity and 
this may complicate attempts to compare optical and X- 
ray lightcurves; for example the optical brightness at a 
fixed mass flow rate, and hence the optical to X-ray flux 
ratio, will depend upon the size of the hot region of the 
disk. 

To perform more useful tests requires UV data. Ideally 
full optical-UV coverage can be used, but the UV alone 
has some power to discriminate. All spectra flatten in the 
UV, although depending on the source and epoch this ef- 
fect may be a modest reduction in slope, or a complete 
turnover (i.e. decreasing F^). Uncertainties in deredden- 
ing the spectra complicate the interpretation of the UV 
slope, but in XTE J1859-I-226 evolution was seen from 
one form to the other, indicating that there are real varia- 
tions. In this case, the evolution was in the opposite sense 
to the optical; the UV spectrum hardened during the out- 
burst decay. These spectral shapes can be thought of as 
defining optical-UV spectral states in BHXRTs, analo- 
gous to X-ray spectral states. The case where the UV 
spectrum continues to rise can be termed the UV-hard 
state, and both A 0620-00 and GS 1124-684 exhibit this 
state. Spectra which peak and fall off in the UV can 
be termed UV-soft, such as those of GRO J0422-h32. In 
this terminology, XTE J1859-I-226 exhibited a transition 
from the UV-soft to UV-hard state as it decayed. 

These spectra can be interpreted in terms of simple 
black body disk models. The full range can be accom- 
modated with two forms of heating, one varying as 
and one as R^'^. The first is readily interpreted as due 
to irradiation by a compact X-ray source and produces 
UV-soft spectra. The second is usually associated with 
viscous heating, but can also arise if the disk is illumi- 
nated by a vertically extended X-ray source. Evolution 
from the first form to the second as the disk cools re- 
sults in exactly the behavior seen in XTE J1859-I-226 - 
the optical spectrum softens, but the UV hardens. None 
of the other sources exhibit such a clear transition, al- 
though GRO J0422-f32 exhibited similar optical spec- 
tral evolution to XTE J1859-I-226 suggested it might 
have shown the same effect if suitable late-epoch data 
were available. A 0620-00 and GS 1124-684 suggest R'^ 
heating, whereas GRO J0422+32 can be interpreted as 
R^^ irradiated heating. It might be significant that both 
GRO J0422-h32 and XTE J185 9-H226 (at least early in 
the outburst) were identified bv lBrocksopp et eH] l|2002jll 
as particularly hard sources; it may be that there is a 
correlation between the optical-UV state (i.e. the form 
of heating of the outer disk) and the X-ray state. Ev- 
idence for this effect is provided by lEsin et all l|2000l^ 
who argue that the transition to the low/hard state near 
the end of the outburst is accompanied by a substantial 
enhancement in the reprocessing efficiency. Obviously a 
larger sample is needed to robustly test these issues. Ide- 
ally we would like to see more examples of transitions in 
the optical-UV spectral state. 

The interpretation of the UV-hard state remains uncer- 



tain. These spectra are well fitted by a viscously heated 
Shakura-Sunyaev disk spectrum. In the case of A 0620- 
00 this interpretation is quite sensible - assuming a lobe- 
filling disk accreting at 30 % of the Eddington rate quite 
naturally produces both the optical fluxes and the spec- 
tral shape without invoking reprocessing as the dominant 
source of optical emission. To explain GS 1124-684 in the 
same way, however, appears to require super-Eddington 
accretion (c.f. ^isra 199^. Irradiation by an extended 
X-ray source offers a way out of this difficulty, and could 
bypass self-shielding of the disk ( Dubus et al. 1999 ). but 
creates its own problems by requiring a very extended 
X-ray emission (or scattering) region. 

Further progress in this area undoubtedly re- 
quires more observations of the quality of those of 
XTE J1859+226. Just as important, however, may be for 
the observations to make contact meaningfully with sim- 
ulations of d isk outbursts. Progress is being made (e.g. 
iLasotal I2001D but simulations have not yet adequately 
accounted for disk warping, or realistic X-ray irradia- 
tion (Dubus et al. 1999); the X-ray emission geometry 
(compact vs. extended), and spectral changes may have 
a significant imp act on the strength of irradiation (e.g. 
lEsin et al.ll2000|) . 
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APPENDIX 

IRRADIATION FROM A CENTRAL SCATTERING CORONA 

We can ex plore different irradiation regime s qualitatively with a simple geometric model of reprocessing, similar to 
that used bv lDumont & Collin-Soviffrin (1990) for active galactic nuclei. Consider a disk element at radius R, height 
H, and local slope dH/dR. The element is illuminated by an isotropically emitting X-ray point source at height Z 
above the center of the disk. The flux F incident upon the surface element depends upon both the distance from the 
X-ray source and the angle of incidence of the X-rays. In general, 

F « ^ (Al) 

where the distance of the X-ray source from the element is d = (— i?, Z — H) in cylindrical coordinates, and the normal 

— 1/2 

to the surface is n = (l + (di//di?)^) {—dH/dR, 1). For a geometrically thin disk (H <g; i?), compact X-ray source 
{Z R), and shallow disk gradients {dH/dR ^1), this result simplifies to: 

F oc ^(dg/di^) + (^-^) (A2) 

The simplest geometry is a flat disk {H = dH/dR — 0). In this case, F oc Z/R^] this is the simple 'lamppost' case 
referred to above. The same would apply if the disk is not flat, but dH/dR is constant, since then H — RdH/dR. 
Somewhat more realistic cases can be approximated with a power-law disk: H = Ho{R/ Rq)^ , where Hq and Rq refer 
to values at the outer disk edge. In this case, 

{I3-1)H + Z 

F oc — (A3) 

/? = 1 is the case described above, and separates the behavior into two regimes. Firstly consider a concave disk, 
P > I. li {P - 1) H <^ Z then the curvature of the disk is negligible and F oc Z/R^ as for a flat disk. If (/? - 1) 7J > Z 
then the vertical elevation of the X-ray source is negligible and we obtain the classical behavior of a concave disk 
illuminated by a central point source, with F oc R~^ in the limiting case of infinitesimal curvature (i.e. /3 — 1 <C 1). 
The dividing line between the two thus depends upon both the height of the disk relative to the height of the X-ray 
source and the amount of curvature of the disk. For an extended X-ray source, radius Z, we might then expect a 
UV-hard spectrum (corresponding to F oc R^^) if the X-ray source has size larger than (/? — 1) H. 

Next consider a convex disk, /? < 1. In this case one never expects heating to vary as R~^. li {1 — f3) H > Z then 
the X-ray source is below the local horizon of the surface element and no irradiation is possible, li {1 — P) H ^ Z 
then again disk curvature is negligible and F oc Z/R^ . In the intermediate cases the radial dependence is even steeper. 
The caveat to this occurs if there is also local scattering material above the disk (i.e. the disk atmosphere) which can 
scatter irradiation over the disk horizon and illuminate regions that would otherwise be obscured. In this case, R~^ 
heating remains possible. 
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